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ABSTRACT 


Some  consequences  of  a  modified  repair  system  for  Phillips'  6^.- 
Appl.  Prob.  IS,  652-659,  Rel .  Engineering,  2,  221-231);iimodel  for  a  two^ 
component  system  are  discussed.  In  the  original  model,  both  components 
are  repaired  whenever  a  revealed  fault  occurs;  in  the  modified  model 

f 

only  faulty  components  are  repaired.  Specifically  (i)  the  distribution 
of  time  from  the  initial  state  up  to  discovery  of  an  unrevealed  fault, 

(ii)  the  e.xpected  proportion  of  time  during  which  there  exists  an 
unrepaired  fault,  and  (iii)  the  distribution  of  number  of  revealed 
faults  up  to  and  including  the  one  which  leads  to  a  discovery  of  an 
unrevealed  fault,  are  obtained.  The  theory  is  illustrated  by  examples, 
based  on  specific  distributions  for  the  times  between  repairs  and  occurrences 
of  the  two  types  of  faults.  A  characterization  of  the  exponential 
distribution  is  indicated. 

Key  Words.  Exponential  distribution,  'Exponential  integral;  Gamma 

distributions.  Geometric  distribution.  Regeneration  point. 
Revealed  faults;  Unrevealed  faults. 
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Introduction 

In  recent  papers,  Phillips  (1979,  1981a, b)  has  investigated  a  system 
consisting  of  two  modules,  faults  in  one  of  which  can  be  revealed  only  on 
special  inspection,  while  the  other  is  monitored  continuously.  A  fault 
developing  in  the  first  module  is  called  'unrevealed'  (U)  (until  a 
special  inspection  is  carried  out),  while  a  fault  in  the  second  module, 
which  it  is  assumed  will  be  detected  immediately,  is  called  'revealed'  (R) . 
It  is  supposed  that  a  special  inspection  for  a  U  fault  is  carried  out  only 
when  a  R  fault  occurs.  Time  for  inspection  and  repair  will  be  neglected 
in  our  analysis,  as  in  our  initial  paper  on  this  subject  (Kotz  and  Johnson 
(1984))  and  in  Phillips'  work. 

Phillips'  model  is  based  on  three  random  variables:  X,  the  time  from 
repair  of  R  to  next  occurrence  of  R  assuming  no  U  present;  Y,  the  time  from 

repair  of  U  to  next  occurrence  of  U,  assuming  no  R  present;  and  Z,  the 

time  from  a  U  fault  to  the  next  R  fault.  X  and  Y  are  assumed  to  be 

independent;  2  is  independent  of  X  but  may  depend  on  Y. 


Phillips  supposes  that  both  modules  are  repaired  whenever  an  R  fault 
occurs  -  whether  there  is  a  U  fault  or  not.  Consequently,  after  each 
repair,  the  system  reverts  to  its  initial  state,  so  there  is  a  regeneration 
point.  Kott  and  Johnson  (1984)  modified  this  model  by  supposing  that  only 
failed  modules  are  repaired.  When  a  U  fault  is  found,  after  an  R 
fault  has  called  for  a  special  inspection,  both  modules  are  repaired 
and  the  system  returns  to  its  initial  state,  (as  it  always  does  in  the 

Phil 1 ips' model ) .  If  a  U  fault  is  not  present  this  does  not  happen,  so 
we  cannot  treat  the  occurrence  of  an  R  fault  alone  as  a  regeneration  point. 

In  either  model,  dependence  of  -  on  Y  makes  allowance  for  the  fact 
that  the  second  module  is  no  longer  new  when  the  U  fault  occurs;  the 
distribution  of  Z  can  also  reflect  the  possible  effects  of  increased  stress 
on  the  second  module  when  the  first  is  not  functioning  properly. 

In  the  present  paper  we  will  (i)  derive  the  distribution  of  time 
(T,  say)  from  initial  state  to  discovery  of  a  U  fault;  (ii)  determine  the 
expected  proportion  of  time  there  is  an  unrepaired  U  fault  in  the  first 
module;  and  (iii)  obtain  the  distribution  of  the  number  of  R  faults  up 
to  and  including  the  one  leading  to  discovery  of  a  U  fault.  As  in  our 
earlier  work,  time  spent  on  repair  will  be  ignored.  (see  Phillips  (I98lb, 
p.  658)  for  the  justification  of  this  assumption).  A  similar  two-component 
model  has  been  analysed  by  Chou  and  Butler  (1985),  assuming  inspection 
to  be  perfect,  in  regard  to  detection  of  defects,  but  possibly  costly 
in  its  effect  on  the  modules.  More  recently,  Murthy  and  Nguyen  (1985) 
also  introduce  cost  considerations  in  a  model  wherein  there  are  no 
unrevealed  faults,  but  occurrence  of  a  fault  in  either  of  the  modules 
may  immediately  cause  a  fault  in  the  other  (with  known  probabilities). 
Repair  is  immediate,  and  only  faulty  modules  are  repaired. 


DISTRIBUTION  OF  TIME  TO  DISCOVERY  OF  A  U  FAULT 


We  will  use  the  same  notation  as  in  Kotz  and  Johnson  (1984)  for 
distribution  functions,  using  the  notation  fjyCw)  to  denote  the  pro¬ 
bability  density  function  (PDF)  of  a  random  variable  W,  and 

00 

\Cw)  =  /  fj^(t)dt 
w 

to  denote  its  survival  function. 

We  also  introduce  f  (x)  to  denote  the  PDF  of  the  n-fold  convolution 

n 

of  f  (x) .  Then 

A 

t 

f^(t)  =  /  S^(y)f,^(y)f,|Y(t|y)dy 

t  y  00 

+  /  /  (I  f  (t)}S  (y-T)f  (y)f7|Y(t-T|y-T)dT  dy 

(t  represents  the  time  of  occurrence  of  the  last  R  fault  preceding  the  R 


expected  number  of  such  events  between  times  a  and  b.) 

It  would  seem  reasonable  to  take 

f_,Y(t:y)=  f*(y+3)/S*(y) 

-00 

where  f*(*)  is  a  PDF  and  S*(y)  =  f*(x)dx  is  the  corresponding  survival 
distribution  function.  This  representation  allows  for  the  age  of  the 
second  module  at  the  time  of  occurrence  of  the  U  fault  in  the  first  module 
in  a  natural  way  (see  e.g.  Kotz  and  Johnson  (1981)).  As  an  illustrative 
example,  take 


•  -v  *7"  *  -  ~'”  •.’"  c* 
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fjjCx) 

a-2  -V0 

=  0  X  e 

o 

A 

<r> 

o 

A 

X 

(1) 

f*(x) 

-2  -x/f 

=  4  X  e 

(x  >  0;  (|)  >  0) 

(23 

(The  expected  values  of  these  distributions  'are  26,  2(J)  respectively , 
usually  we  would  have  6.3  For  th_e  present  we  do  not  specify  fyCy) 


From  (1),  f^Cx3  =  /  (2n-l)! 


-1  ,  -x/i) 


S  (x3  =  (l+e"^x}e‘^''®  ;  s*(x)  =  (!+(?'  x}e 


a 


Also 


I  f  (x)  =  0'^  I  {(X9'^)^"'^  /  (2n-l)!} 

’  n=l  ”  n=l 


Hence, 


f_(t3  =  /  (l+e'^y)e’^'^®f  (y)c:)'^  t  ^dy 

T  n  ^ 


0 

0  0" 


.  f'  Cly)  •  exp{-(S-'-<t-')y)dy 

0  i*y  y 

-  0  0 


,-l, 


.  d-*-9  Cy-T)  (t-T)fY(y)exp{- (0  ^-4)  ^)y+9  ^T}dT  dy  (3) 
1+?  (y-x) 

The  expression  for  f-i-(t3  looks  formidable,  but  can  be  simplified  for 
some  particular  choices  of  f^Cy)- 


‘-y 

I  tb.\ 


PROPORTION  OF  TIME  THERE  IS  AN  UNDISCOVERED  U  FAULT 


At  the  time  of  the  first  R  fault  after  occurrence  of  a  U  fault, 
both  modules  are  repaired,  and  the  system  starts  again  from  its  initial 
state.  So  the  log-run  proportion  of  time  there  is  an  undiscovered  U 
fault  is  the  same  as  the  expected  proportion  of  time  this  is  so  in  a  singl 

cycle  to  discovery  of  the  first  U  fault.  This  is  E[Z]/E[T],  where  T 
denotes  time  to  first  discovery  of  a  U  fault.  (Note  that  it  is  not 
E[Z/T],  which  would  give  undue  weight  to  the  proportions  for  small  values 
of  T.) 

We  can,  in  principle,  evaluate  E[T]  from  the  PDF,  f.j.(t),  as 
00 

E[T]  =  /  t  f.^(tjdt  . 

0 

30 

Since  Z  =  T-Y,  E[Z]  =  E[T]  -  E[Y]  and,  of  course  E[Y]  =  /  y  fY(y)dy 

0 

In  general,  however,  evaluation  of  E[T]  is  difficult  because  of  the 
complexity  of  f.p(t)  (see  (3)).  Considerable  simplification  results  if  Y 
has  an  exponential  distribution 

fyCyJ  (y>0;  ■^>  0)  (4) 

whereby  we  can  avoid  using  f.j,(t3  explicitly.  Because  of  the  lack-of- 
memory  property  of  the  exponential  distribution,  the  system  reverts  to 
its  initial  state  after  any  repair,  and  in  particular  after  the  first 
repair  of  an  R  fault,  whether  or  not  there  is  a  U  fault  to  be 
repaired  at  the  same  time.  So  we  need  only  evaluate  £[Z']/E[T']  where  T^ 
is  the  time  to  first  occurrence  of  an  R  fault,  and  Z'  is  the  time  (which 
may  be  zero)  there  is  a  U  fault  before  any  R'  fault  occurs.  Wic 
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-  L  00 

=  /  /  S_^(y3f^[y)  {f*(t)/S*(y}}  t  dy  dt  +  /  t  SyCtjf^ (t)dt  . 


The  first  and  second  terms  on  the  rtgfit-hand  side  correspond  to 
the  R  fault  being  preceded,  or  not  (respectively)  by  a  U  fault.  Also 
from  (5) 


00  00  oo 

e[t1  =  e[^  +  /  yf^Cy)  {S^(y)/S*(y) }/  f*(t)dt  +  /_  t  SYCt)f^Ct)dt 


=  E[:]  +  /  yifyCyJS^Cy)  +  s  (y)f,^(y)}dy  (6; 

0 

Inserting  the  PPF's  from  (1),  (2)  and  (4),  and  putting  0^+u)^= 
we  have  from  (3) 


E  [^]  =  /  j  (1+  9  ®  *  4>  ^  t  e  (1+<|)  ^y)  e  ^  (t-y)dy  dt 


;  j 

0  0 


=  (-oO")'^  /  exp{-(3~^-Q'^)y}  J  t(t-y)e'^^^dt  dy 

0  1+:  y  y 

00 

=  ^  f  - — ^  •  [2.i^{l+(l)  V  *  iCb  ■  4>^y(l  +  i+>  V)]e  ^'^^dy 

^0  l+:'\v 

-1  1  +  9‘^v  -1  ,  -y/3, 

=  .o  :  j  - (2+0  y}e  dy 

0  1+:  y 


-1  -y/3, 

,  0  '  (l  +  r  y  +  - j^)e  dy 

0  l+o'  y 


{- +  +  b(i-5e'^)Ej(0B'^)e^/-}  (7) 


.  « ^  .  VVX 

where  E^(wi  =  J  x  e  dx  is  the  exponential  integral  (see  Appendix  I). 

1 

00  OO 

Also,  from  (6),  E[T^  =  E[2^  +  J  yy  ^e  ^^"“(l+S  V)®  ^^"dy  +  J  te  “t  e  ^'^^dt 

0  0 


00 

=  E[:^  +  /  {oj  ^y(l+e'\v}  +  9''‘y‘'}e'^^®dy 


=  E[:^  *  v‘^(6"+2e‘^g'')  +  29*“  S'"  =  E[:^  +  3 


/iius  tiie  rati„  pk)',--  is  the  same  as 

^  E(T) 


The  ratio  E[Z]/E[T]  depends  only  on  the  ratios  The  following 

table  gives  a  few  illustrative  values.  It  should  be  remembered  that 
the  expected  values  of  the  three  distributions  (1),  (2)  and  (4)  are 
29,  2b  and  lo  respectively  (not  9,0  and  u)) . 


Table  1.  E[Z]/E[T]  =  Proportion  of  Time  There  is  an  Undiscovered  U  Fault 

a.  /  '^ 


^/e  '  0.4  0.6  O.S  1.0 


0.5 

1 

0.546 

0.654 

0.723 

0.778 

1 

0.550 

0.440 

0.522 

0.625 

2  1 

0.187 

0.270 

0.340 

0.444 

4 

1 

0.095 

0.141 

0.188 

0.280 

8  i 

0.042 

0.066 

0.088 

0.160 

As  is  to  be  expected,  the  longer  the  expected  life  time  of  the  first 
module  (^)  the  less  proportion  of  time  there  is  an  undiscovered  U  fault. 
The  greater  the  expected  value  of  Z,  given  Y  (i.e.  the  greater  i)  the 
greater  the  proportion  of  time  there  is  an  undiscovered  U  fault.  The 
rate  of  increase  is  more  pronounced  for  larger  values  of  jj/ 9  . 


The  model  considered  by  Philips  corresponds  to  our  model  with  the 


addition  of  the  assumption  that  f,iY(z|y)  =  f.,(zj,  independent  of  y. 


DISTRIBUTION  OF  NUMBER  OF  .R  FAULTS  UP  TO  A.N’D  INCLUDING  DISCOVERY  OF  A  U  FAULT 
Denoting  the  number  of  R  faults  up  to  and  including  discovery  of  a  U 


fault  by  N,  we  have 
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Pr[N=„J  =  f,,Cy){S^(y)  -  S__j(y);.dy  -  E^[S_^(y)  -  S^_j(y))  (8) 

00  y 

where  S  (y)  =  /  f  CyJdy  =  1  -  [  f  Cy)dy 

y  n  0  " 

Integrating  by  parts,  we  obtain  the  equivalent  formula 

OO 

Pr[N=n]  =  /  S  (y){f  (y)  .  f  (y) }dy  (8)' 

Q  ^  ii  1  Jj 

In  particular,  in  the  exponential  case,  with 
fyCy)  =  (y>0;  >  0) 

then  because  of  the  lack-of-memory  property  of  the  exponential  distribution, 
the  distribution  of  N’  is  geometric  with 

Pr[N--n]  =  Pd-Pj""’^  (n=l,2,...)  (9) 


where  P  =  Pr[Y''X].  Subject  to  the  condition  that  the  distribution  of 

Yr>0)  is  absolutely  continuous,  the  converse  is  true.  That  is,  if  a 

distribution  of  form  (9)  for  N  holds,  whatever  be  the  density  function 

f^(x),  then  Y  must  have  an  exponential  distribution  (cf.  Grosswald  et  al 

(1980)).  A  proof  is  presented  in  Appendix  II. 

If  X  has  the  distribution  (1)  then 


1  -  P 


00 

i 

0 


dx 


so 


Pr[N=n]  -  (1-- 


r)(e' 


.2,n-l 


(n=l,2,...) 


In  terms  of  u,  6  “  d"  =  9  ~ (S  ^+'0)'^)  “  =  uj“(6+w)'‘' 


Formula  (3)  is  quite  general.  For  example,  if 


y  (y>0;  ^  >  0) 


(10) 
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Using,  in  CD  and  (lOJ,  general  Erlang  distributions,  we  can  generate 
systematically,  via  (S},  a  family  of  discrete  analogues  of  these  distributions, 
of  which  (11)  is  a  particular  case. 
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APPENDIX  I 


/  (l+at)(l+bt)  =  b  ^  J  (l+ab* ^u)  C  1+u)  ” 

0  0 

,-l  c/b  f  I-  i -I  T  -I  -cv/b, 

=  b  e  j  {1+ab  (v-l)}v  e  dv 

1 

(with  v=l+u) 

u-1  c/b  r°°  ,  .  -1  ,,  ,-1-1  -1,  -cv/b 

-b  e'  j  {ab  +(l-ab  jv  }e  dv 

1 

=  b'^{ac'^+(l-ab'hE^(ab'^)e‘^'^^}  . 

OO 

where  E  (oj)  =  -EiC-uj)  =  /  x  e  ^^dx  is  the  exponential  integral  (tabulated  ^ 
for  example,  in  U.S.  National  Bureau  of  Standards  (1941)), 

Putting  a=5  \  b=i  ^  and  c=3  ^  we  obtain 

I  ^  e'^^^dt  =  y{3e'^  +  (l-4»0'^)E  ((?3"^)e^'^'^} 

This  formula  is  used  in  equation  (7). 


-12- 


APPENDIX  II 

We  will  show  that  if  the  distribution  of  N  is  of  form  (9),  whatever 
the  density  function  f^Cx),  and  Y  is  absolutely  continuous,  then  Y  must 
have  an  exponential  distribution. 

Take  the  special  case  fj^(x)  =  e-^  (x  >  0)  whencef^^Cx)  =  x"‘^e'^/ Cn-1 ) ! 

and  S^(y)-S^_^(y)  =  y"'^e‘^/ Cn-1) ! .  C1-) 

If  g^(y)  and  g^Cy)  are  two  different  density  functions  of  Y  giving  the 
same  formula  for  Pr[N=n]  then  from  (8)  and  (12) 

{e‘^g^(y) ly^'^dy  =  /^  {e  ^g2(y)}y^  ^dy  (n=l,2,,..) 


This  means  that  the  two  density  functions  c^.  e  ^gj  (y)  (0<y)  with 

c.  =  {j“  e'^g^ Cy)dy}'^  (j=l,2),  have  the  same  moments  of  positive  integer  order. 

Since  the  exponential  density  (e  is  determined  by  its  moments,  and 
g^(;y3,g^(y)  are  density  functions  the  density  functions  c^e‘^g^(y)  are 
determined  by  their  moments,  whence 

c^e'-^'g^Cy)  2  c^e'^g^Cy) 


so  g^Cy)  =  g2Cy)  (and  Ci=C2)- 

Since  gCy)  =  f^Cy)  =  (0<co;0<y)  does  satisfy  (9),  it  follows 

that  g^(y)  (g-,(y))  must  be  of  this  exponential  form. 

The  result  clearly  holds  under  the  broader  conditions  that  (8)  holds  for 

some  S  (y)  -  S  Ay)  =  y"'^h(y)  where  h(y)/(/^  h(y)dy}  is  a  density  function 
n  n- 1  V 

determined  by  its  moments. 
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